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ABSTRACT 



Structural brain abnormalities identified at near-term age have been recognized as potential predictors of 
neurodevelopment in children born preterm. The aim of this study was to examine the relationship between neo- 
natal physiological risk factors and early brain structure in very-low-birth-weight (VLBW) preterm infants using 
structural MRI and diffusion tensor imaging (DTI) at near-term age. 

Structural brain MRI, diffusion-weighted scans, and neonatal physiological risk factors were analyzed in a cross- 
sectional sample of 102 VLBW preterm infants (BW < 1500 g, gestational age (GA) < 32 weeks), who were ad- 
mitted to the Lucile Packard Children's Hospital, Stanford NICU and recruited to participate prior to routine near- 
term brain MRI conducted at 36.6 ± 1.8 weeks postmenstrual age (PMA) from 2010 to 2011; 66/102 also 
underwent a diffusion-weighted scan. Brain abnormalities were assessed qualitatively on structural MRI, and 
white matter (WM) microstructure was analyzed quantitatively on DTI in six subcortical regions defined by 
DiffeoMap neonatal brain atlas. Specific regions of interest included the genu and splenium of the corpus 
callosum, anterior and posterior limbs of the internal capsule, the thalamus, and the globus pallidus. Regional 
fractional anisotropy (FA) and mean diffusivity (MD) were calculated using DTI data and examined in relation 
to neonatal physiological risk factors including gestational age (GA), bronchopulmonary dysplasia (BPD), necro- 
tizing enterocolitis (NEC), retinopathy of prematurity (ROP).and sepsis, as well as serum levels of C-reactive pro- 
tein (CRP), glucose, albumin, and total bilirubin. 

Brain abnormalities were observed on structural MRI in 38/102 infants including 35% of females and 40% of 
males. Infants with brain abnormalities observed on MRI had higher incidence of BPD (42% vs. 25%) and sepsis 
(21% vs. 6%) and higher mean and peak serum CRP levels, respectively, (0.64 vs. 0.34 mg/dL, p = .008; 1.57 vs. 
0.67 mg/dL, p= .006) compared to those without. The number of signal abnormalities observed on structural 
MRI correlated to mean and peak CRP (rho = .316, p = .002; rho = .318, p= .002). The number of signal abnor- 
malities observed on MRI correlated with thalamus MD (left: r= .382, p= .002; right: r= .400, p= .001 ), control- 
ling for PMA-at-scan. Thalamus WM microstructure demonstrated the strongest associations with neonatal risk 
factors. Higher thalamus MD on the left and right, respectively, was associated with lower GA (r = —.322, p = 
.009; r= -.381, p= .002), lower mean albumin (r= -.276, p= .029; r= -.385, p= .002), and lower mean bil- 
irubin (r = -.293, p= .020; r= -.337 p= .007). 

Results suggest that at near-term age, thalamus WM microstructure may be particularly vulnerable to certain 
neonatal risk factors. Interactions between albumin, bilirubin, phototherapy, and brain development warrant fur- 
ther investigation. Identification of physiological risk factors associated with selective vulnerability of certain 
brain regions at near-term age may clarify the etiology of neurodevelopmental impairment and inform neuro- 
protective treatment for VLBW preterm infants. 

© 2014 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY-NC-SA license 

(http://creativecommons.Org/licenses/by-nc-sa/3.0/). 



Abbreviations: VLBW, very-low-birth-weight; GA, gestational age; PMA, post-menstrual age; DTI, diffusion tensor imaging; FA, fractional anisotropy; MD, mean diffusivity; CC, corpus 
callosum; IC, internal capsule; ALIC, anterior limb of the internal capsule; PL1C, posterior limb of the internal capsule; GloP, globus pallidus. 
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1. Introduction 

Children born preterm with very-low-birth-weight (VLBW) are at 
higher risk for neurodevelopmental impairments, such as cerebral 
palsy, developmental delay, and possibly autism spectrum disorder 
(Woodward et al., 2012; Movsas et al., 2013). Brain abnormalities iden- 
tified in the neonate may result in higher risk of neurologic impairment 
and later functional deficits (Woodward et al., 2012). Rapidly develop- 
ing early brain regions may be selectively vulnerable to injury that 
may underlie adverse outcome (Johnston and Hoon, 2000; Miller and 
Ferriero, 2009; Liu et al., 2013). Previous studies have demonstrated 
widespread differences in white matter (WM) development and re- 
gional cortical volumes in infants born VLBW compared to age- 
matched full-term infants (Hiippi et al., 1998; Inder et al., 2005; 
Thompson et al., 2007; Lee et al., 2013; Thompson et al., 2013). These 
structural differences may be biomarkers for early brain injury and ad- 
verse outcome. 

Near-term neuroimaging data can inform early prognosis for pre- 
term infants and has potential to guide early intervention at a time of 
optimal neuroplasticity. Specifically, structural MRI can inform about 
gross brain injury, which has been shown to correlate with outcome 
(Miller and Ferriero, 2009; Liu et al., 2013); however, diffusion tensor 
imaging (DTI) of early WM regions may be a more sensitive biomarker 
of later neurodevelopment than structural MRI (Arzoumanian et al., 
2003; Rose et al., 2007; Rose et al., 2009; De Brume et al., 2013). DTI 
quantifies the diffusion properties of tissues and is used to investigate 
WM microstructure. As the brain develops, water content decreases, ex- 
tracellular spaces diminish, and cellular microstructures become more 
complex and organized; these changes can be quantified with DTI 
(Dubois et al., 2008; Lee et al., 2013; Dubois et al., 2013; Oishi et al., 
2013; Nossin-Manor et al., 2013). 

Scalars obtained from DTI can assess brain development and matu- 
ration and include fractional anisotropy (FA), mean diffusivity (MD), 
axial diffusivity (AD), and radial diffusivity (RD) (Pierpaoli et al., 
1996). FA reflects the degree of diffusion anisotropy within a voxel, 
and is determined by fiber diameter and density, myelination, extracel- 
lular diffusion, inter-axonal spacing, and intra-voxel fiber-tract coher- 
ence. MD is a calculation of average diffusion along the three main 
axes, relative to the primary direction of diffusion, or AD. AD is thought 
to reflect fiber coherence and structure of axonal membranes. RD, the 
mean of diffusivities perpendicular to the primary axis, is thought to 
represent the degree of myelination (Song et al., 2002; Chen et al., 
201 1 ). In neonatal WM, FA has been found to increase with age, while 
MD, AD, and RD decrease with age, likely due to increased fiber organi- 
zation, axonal coherence, and preliminary myelination (De Brui'ne et al., 
2013; Rose et al., 2014). Analyses of DTI scans obtained during the neo- 
natal period may serve as potential predictors of neurodevelopment in 
children born preterm (Hiippi and Dubois, 2006; Rose et al., 2009; 
Thompson et al., 2013). 

Previously, we reported DTI findings on regional brain microstruc- 
ture at near-term age in 45 VLBW preterm infants with no evidence of 
gross brain abnormalities, a subpopulation of infants reported in the 
current study (Rose et al., 2014). WM microstructure was analyzed in 
19 subcortical regions. Results indicated a general pattern of central to 
peripheral regional WM development and of posterior to anterior WM 
development within regions of the corona radiata, corpus callosum 
(CC) and internal capsule (IC). Projection fibers of the PLIC demonstrat- 
ed the highest degree of WM maturation on DTI (Vassar et al., 2013; 
Rose et al., 2014). The ALIC was less mature but appeared to be one of 
the most rapidly developing WM regions at near-term age, consistent 
with early stages of myelination. White matter within the thalamus 
and globus pallidus (GP) demonstrated higher FA and lower MD values 
and had larger, less-variable volumes selected, than the putamen and 
caudate (Rose et al., 2014), suggesting that these regions may be more 
reliably assessed at near-term age. These findings, as well as other re- 
cent studies (Dubois et al., 2013), suggest that the developmental status 



of WM regions is an important consideration in the interpretation of 
neonatal neuroimaging studies, and informed the selection of regions 
analyzed in the current study. 

Rapidly developing brain regions in the neonate may be selec- 
tively vulnerable to injury due to morbidities and clinical factors 
common during the neonatal period. Consequently, abnormal WM 
development in these regions may be useful for identifying infants 
with higher risk for adverse outcomes. Known risk factors for poor 
neurodevelopment following preterm birth include sepsis (Stoll 
et al., 2002; Schlapbach et al., 2011; Hentges et al., 2013; Mitha 
et al, 2013), bronchopulmonary dysplasia (Schmidt et al., 2003; 
Vohr et al., 2005; Karagianni et al., 2011; Schlapbach et al., 2012), 
necrotizing enterocolitis (NEC) (Hintz et al., 2005; van Vliet et al., 
2013), and male sex (Wood et al., 2005; Vohr et al., 2005; Rose 
et al, 2009; Beaino et al., 2010; Morsing et al., 2011). The aim of 
this study was to examine the relationship between neonatal physi- 
ological risk factors and early brain microstructure in VLBW preterm 
infants at near-term age, based on qualitative analysis of structural 
MRI and on semi-automated, atlas-based DTI analysis. 

2. Methods 

2.1. Subjects 

Near-term scans were obtained from 102 VLBW neonates, 
representing 76% of eligible participants treated at Lucile Packard 
Children's Hospital (LPCH) between 1/1/10 and 12/31/11. Parents 
were approached prior to scheduled MRI and consent was obtained 
for this IRB-approved study. Infants with GA-at-birth <32 weeks, birth 
weight <1500 g, and no evidence of genetic disorders or congenital 
brain abnormalities were included. Sixty-six of the 1 02 infants had suc- 
cessful DTI scans, collected at the end of routine near-term MRI scans, 
prior to discharge. 

2.2. Neonatal risk factor definition and data collection 

Demographic data and neonatal physiological risk factors were 
assessed in 102 VLBW preterm infants as recorded in the medical re- 
cord. Demographics included gestational age (GA), birth weight (BW), 
sex, Apgar score at 5 minutes postnatal (Apgar-at-5 min), multiple 
births, maternal age, pregnancy-induced hypertension (PIH), and pre- 
mature rupture of membranes. Physiological risk factors included the 
presence of bronchopulmonary dysplasia (BPD) defined as having a his- 
tory of respiratory distress syndrome and treated with oxygen at 
36 weeks PMA, sepsis confirmed with positive blood culture (Sood 
et al., 2012), necrotizing enterocolitis (NEC), and retinopathy of prema- 
turity (ROP) at any time during hospital stay. Serum levels of mean and 
peak C-reactive protein (CRP), as well as mean glucose, albumin, and 
total bilirubin levels during the first 2-weeks of life were recorded. Lab- 
oratory values were obtained through a medical record search conduct- 
ed with Stanford Center for Clinical Informatics and STRIDE (Stanford 
Translational Research Integrated Database Environment) (Lowe et al., 
2009). 

Neonatal demographic and physiological risk factors were selected 
for analysis based on their association with poor neurodevelopmental 
outcome in previous studies, as detailed in the Introduction. As noted 
in Table 1, serum levels of CRP, glucose, albumin, and total 
(T) bilirubin were sampled multiple times over the first two postnatal 
weeks. Other risk factors were identified based on physician report in 
the medical record. 

2.3. MRI data acquisition 

Brain MRI scans were performed on 3T MRI (GE-Discovery MR750, 
GE 8-Channel HD head coil) at Lucile Packard Children's Hospital at 
Stanford University. The MRI included Tl-weighted, T2-weighted, and 
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Table 1 

Demographic and physiological risk factors. 





All Npnnat"pc 

n = 102 


Nnnnal MRI 
iNuiiiuu lvirvi 

n = 64 


Ahnnmnal MRI 

/ MJl IU 1 1 1 l<l 1 Ivllvl 

n = 38 


MpAnafpc witTi nTl 

lNCUlldLC3 VV1L11 \J 1 1 

n = 66 


Males n (% with abnormal MRI) 


42 


25 


17 (40%) 


25 


Females n (%with abnormal MRI) 


60 


39 


21 (35%) 


41 


Gestational age (weeks) mean (SD) 


28.7 (2.4) 


28.9 (2.3) 


28.3 (2.6) 


28.9 (2.3) 


Birth-weight (g) mean (SD) 


108 (279) 


1102 (276) 


1162 (285) 


109 (266) 


PMA-at-scan (weeks) mean (SD) 


36.6 (1.8) 


36.5 (1.6) 


36.8 (2.0) 


36.5(1.3) 


Maternal age (years) mean (SD) 


31.6 (6.0) 


32.2 (5.4) 


30.6 (6.8) 


31.9(6.1) 


Multiple births mean (SD) 


1.75 (.95) 


1.8 (.88) 


1.7(1.1) 


1.8 (.92) 


Apgar at-5-min mean (SD) 


7.4(1.8) 


7.7 (1.6) 


6.8 (2.2) 


7.4 (1.7) 


Days on ventilation mean (SD) 


11.1 (18.4) 


7.7 (13.9) 


16.7 (23.2) 


8.7 (14.4) 


R rnnrhniu ilnihni rw H wcr^l 3d 3 n f 
Dl UllLllUpUllllUlldl y UybpidMd, 11 \/o) 


32 (31%) 


1 fi f?5%1 


16 (42%) 


18 (27%) 


Necrotizing enterocolitis, n (%) 


14(14%) 


8 (13%) 


6 (16%) 


7 (11%) 


Retinopathy of prematurity, n (%) 


72 (71%) 


40 (63%) 


32 (84%) 


45 (68%) 


Sepsis a , n (%) 


12 (12%) 


4 (6%) 


8 (21%) 


59 (9%) 


Mean CRP b , mg/dLmean (SD) 


.45 (.58) 


.34 (.45) 


.64 (.72/ 


.38 (.49) 






p = .008 






Peak CRP b , mg/dL mean (SD) 


1.00 (1.49) 


.67 (.98) 


1.57 (1.99)' 


.83 (1.30) 






p = .006 






Mean glucose 1 , mg/dL mean (SD) 


103.1 (20.9) 


101.8 (19.3) 


105.2 (23.7) 


103.0 (18.3) 


Mean albumin' 1 , mg/dL, mean (SD) 


2.4 (.3) 


p = .488 
2.4 (0.3) 


2.4 (0.4) 


2.4 (.3) 






p = .185 






Mean T bilirubin 6 mg/dL, mean (SD) 


5.4(1.3) 


5.4 (1.3) 


5.4 (1.4) 


5.3 (1.3) 






p = .802 







a Confirmed with a positive blood culture. 

b CRP sampled an average of 3.72 (1.9) times in the first 2-weeks. 

c Glucose sampled an average of 13.6 (5.8) times in the first 2-weeks. 

4 Albumin sampled 9.7 (6.1) times in the first 2-weeks. 

e Total bilirubin sampled 13.6 (2.9) times in the first 2-weeks. 

f Significantdifference between infants with normal versus abnormal MRI assessed with Mann-Whitney U using two-tailed significance of p < .01 for serum CRP, glucose, albumin, and 
bilirubin. 



diffusion-weighted imaging (DWI) and diffusion tensor imaging (DTI) 
scans. A gradient echo 3-plane localizer was used, and an asset calibra- 
tion was proscribed to utilize parallel imaging. From the 3-plane 
localizer, Tl images were generated with the following parameters: 
TE = 91.0, TR = 2200, FOV = 20 cm, matrix size = 320 x 224, slice 
thickness 3.0 x 0.5 mm spacing, NEX = 1. T2, DWI, and DTI axial 
scans were proscribed using a single acquisition, full phase field of 
view (FOV). The Fast Recovery, Fast Spin Echo T2 scan imaging param- 
eters were as follows: TE = 85 ms, TR = 2500, FOV = 20 cm, matrix 
size = 384 x 224; slice thickness = 4.0 mm x 0.0 mm spacing. Axial 
T2 FLAIR parameters were as follows: TE = 140, TR = 9500, FOV = 
20 cm, slice thickness = 4.0 mm x 0.0 mm, inversion time = 2300, 
fluid attenuated inversion recovery matrix = 384 x 224. Axial DWI pa- 
rameters were as follows: TE = 88.8, TR = 10,000, FOV = 20 cm, slice 
thickness = 4.0 mm x 0.0 mm spacing, matrix = 128 x 128. The DTI 
scan was based on a diffusion-weighted, single-shot, spin-echo, echo- 
planar imaging sequence with a slice thickness of 3.0 mm, a matrix 
size of 128 x 128, a 90° flip angle, FOV = 20 cm, TE = 88.8 ms, TR = 
8000, (b = 1000 s/mm 2 ). Diffusion was measured along 25 directions, 
with three B0 images. Two repetitions were obtained from 64/66 sub- 
jects. MRI scans were performed as the routine near-term neuroimaging 
for preterm infants and the diffusion-weighted sequences were collect- 
ed at the end of these ~25-minute MRI scans. Infants were swaddled and 
fed and typically remained asleep for the duration of the scan. Sedation 
was not utilized as part of the research protocol and typically was not 
utilized for routine near-term MRI, although it may have been pre- 
scribed according to clinical needs in some cases. 

2.4. Radiological assessment 

Brain abnormalities were assessed on structural MRI for degree of 
abnormality on a scale that defines categories of 1-5, as described in 
Table 2, and for the number of signal abnormalities observed on Tl 



Table 2 

Categories used to classify structural MRI findings. 



Category 


Radiological findings on MRI 


CI 


No abnormality. 


C2 


Minimal abnormality, such as subependymal hemorrhage or 




mineralization 




(e.g. sequelae of germinal matrix hemorrhage) with or without mild 




ventriculomegaly, with or without mild IVH. 


C3 


Moderate to severe ventriculomegaly, with or without IVH 


C4 


Parenchymal abnormality, including evidence of gray and/or white matter 




injury due to hemorrhage or ischemia, with or without IVH. 


C5 


Congenital malformation (excluded from further analyses) 



and T2 scans in cerebral, intraventricular, internal capsule, basal ganglia, 
thalamus, and cerebellar regions. Radiological evaluation was per- 
formed by an experienced pediatric neuroradiologist, X.S., who was 
masked to all other participant data. Brain abnormalities identified in 
the research assessment were consistent with findings of brain abnor- 
malities reported by the clinical neuroradiologist in the patient's medi- 
cal record. 



2.5. DTI processing 

In a subpopulation of 66 infants, semi-automated, atlas-based DTI 
was analyzed. Diffusion-weighted images were pre-processed using 
DTI-Studio based on parameters outlined in Jiang et al. (2006) and 
Oishi et al. (2011). For infants with two repetitions (n = 64/66), the 
best repetition was selected to eliminate images with artifacts or evi- 
dence of motion. If no full repetition was usable, a composite repetition 
was generated based on best image slices (provided the two repetitions 
were well-registered). Automated Image Registration was performed 
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Fig. 1. Brain regions selected in a representative infant DTI scan using JHU neonatal atlas, 
axial and corona views. 



using an affine transformation to further correct for eddy current distor- 
tions and minor head motion. 

Each brain image was skull-stripped with ROl Editor using BO and 
trace images to distinguish brain matter from intracranial CSF and 
skull, and manually rotated to align with the JHU neonatal template 
image. DTI images were processed with DiffeoMap (http://www. 
mristudio.org ) using FA and trace (the sum of the 3 orthogonal diffusiv- 
ities per voxel) maps to perform a large deformation diffeomorphic 
metric mapping (LDDMM) transformation (Oishi et al., 2011). Each 
brain was normalized to map onto the neonatal atlas and segmented 
into 126 regions. Based on the trace map, regions >0.006 mm 2 /s were 
considered CSF and were excluded from the segmented regions. A 
threshold of FA >0.15 was then applied to segment developing WM 
tracts (Vassar et al., 2013). Regions selected within the FA >0.15 thresh- 
old were then mapped onto the FA and trace maps, and regional values 
were obtained. MD values were defined as an average of the three direc- 
tions of diffusion calculated from the trace image. 

We examined subcortical WM defined by the JHU panellation atlas 
that was well-visualized with near-term MR1 and DTI, and known to 
mediate motor or cognitive function, as illustrated in Fig. 1 . Specifically, 
we measured FA and MD in commissural fibers of the genu and 
splenium of the ccs, projection fibers of anterior and posterior limbs of 
the internal capsule (ALIC and PL1C respectively) and in the WM of sub- 
cortical nuclei of the thalamus and globus pallidus (GloP). 

These regions were selected for analysis based on prior work that 
demonstrated that these regions were reliably detected by semi- 
automated, atlas based analysis of DTI at near-term age (Rose et al., 
2014) based on percent volume above the FA >0.15 threshold, relative 
to original volumes identified by the neonatal atlas, as applied to the 
45 infants with no evidence of abnormalities on MRI and scanned 
prior to 40 weeks PMA (Vassar et al., 2013). 



2.6. Statistical analysis 

Neonatal demographic and physiological risk factors of GA, CRP, glu- 
cose, albumin, and bilirubin were compared in infants with and without 
brain abnormalities observed on MRI using nonparametric Mann- Whit- 
ney U, with a significance level of p< .01 as a correction for multiple 
comparisons (Table 1 ). The Spearman's correlation was used to assess 
associations between neonatal risk factors and brain MRI score and 
the number of signal abnormalities observed on MRI with a significance 
level of p< .008, to correct for multiple comparisons (Table 3). Partial 
correlations were reported between MRI score and number of signal ab- 
normalities observed on MRI and MD values, after controlling for PMA- 
at-scan with a significance level of p< .01 as a correction for multiple 
comparisons. Partial correlations were reported between neonatal risk 
factors and MD values in the CC, ALIC, PLIC, thalamus, and GloP, after 
controlling for PMA-at-scan with significance level of p < .01 —corrected 
for multiple comparisons (Table 4). Differences between male and fe- 
male MD values, corrected for PMA-at-scan, were compared using the 
nonparametric Mann-Whitney test, with a significance level of 
p < .05. Spearman's correlation was used to assess associations among 
risk factors with continuous values with a significance level of p< .01 
(Table 5). The variance inflation factors (VIF) were estimated as a diag- 
nostic for co-linearity between risk factors. Statistical analysis was per- 
formed using SPSS Statistics software Version 21 (IBM). 

3. Results 

Demographic characteristics and physiological risk factors of all par- 
ticipants who had MRI (n = 102) and the subset population that also 
had DTI (n= 66) are reported in Table 1. Demographic data of the 
total population was similar to that of the sub-population who had 
both MRI and DTI. Demographic and physiological risk-factor data for 
participants with and without brain abnormalities observed on MRI 
are reported in Table 1 . Brain abnormalities were identified on MRI in 
38 out of 102 participants (37%), including 21/60 (35%) females and 
17/42 (41%) males. Gestational age and BW were similar in those with 
and without brain abnormalities on MRI. Infants with brain abnormali- 
ties observed on MRI had higher incidence of BPD (42% vs. 25%) and 
sepsis (21% vs. 6%) and higher mean and peak serum CRP levels (0.64 
vs. 0.34 mg/dL, p= .008; 1.57 vs. 0.67 mg/dL, p = .006) compared to 
those without brain abnormalities observed on MRI. 

The correlations between neonatal risk factors are shown in Table 5. 
Of note, a strong association was found between mean bilirubin and BW 
(rho = .789, p = .000) and between mean albumin and mean bilirubin 
levels (rho = .604,p = .000). For risk factors with continuous values we 
estimated the variance inflation factors (VIF) as a diagnostic for collin- 
earity and found that mean and peak values of CRP were collinear 
(VIF < 9.9). Otherwise, among the mean values of risk factors we 
found no other collinearities, indicating that these mean values provide 
reasonably independent information for predictive models. The VIF 
ranged from 1.0 to 3.2, and fell below 10, the conservative criterion for 
indicating presence of collinearity. 

Table 3 lists the mean and range of MRI score and the number of sig- 
nal abnormalities observed on MRI, as well as correlations with physio- 
logical risk factors. The higher number of signal abnormalities was 
correlated to lower scores for Apgar-at-5-minutes and higher mean 
and peak CRP (r = .316, p= .002; r= .318, p= .002). 

DTI analysis was performed on the selected brain regions shown in 
Fig. 1 as defined by the JHU neonatal atlas. Detection of brain abnormal- 
ities observed qualitatively on structural MRI compared to DTI findings 
indicated that the number of signal abnormalities observed on MRI was 
associated with thalamus MD (left: r= .382, p = .002; right: r= .400, 
p = .001 ), with trends for genu, splenium, ALIC, and PLIC MD. No asso- 
ciations were found between signal abnormalities observed on MRI and 
GloP MD. 
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Table 3 

Correlations between neonatal risk factors and brain MR1 score and number of signal abnormalities on MR1 at near-term age. 



1NCU1 IcLlcLl 1 1 jl\ IcLL LU1 :> 




R3Ilg6 — 1—4 

iVlCdll ^JUJ ^ lt\J) 


it nf cional alinnrmalitipc 

Tr Ul M^llm dUl 1U1 11 kill LICo 

Piim-a — n 1 ^ 
Kdllgc — U— 1 J 

ypin fern — 91 f3 f]) 

IvlCall ^3L/ J — Z., 1 ^ J-L* J 


Gestational age at birth 


Rho 


-.140 


-.197 




Sig. 


p = .163 


p = .048 


Apgar at 5 min 


Rho 


-.236 


-.303 




Sig. 


p = .031 


p= .005 a 


C-reactive protein 


Rho 


.220 


.316 




Sig. 


p = .031 


p= .002 a 


Peak C-reactive protein 


Rho 


.211 


.318 




Sig. 


p = .039 


p= .002 a 


Glucose 


Rho 


-.001 


-.028 




Sig. 


p = .994 


p = .787 


Albumin 


Rho 


- .170 


-.193 




Sig. 


p = .095 


p = .058 


Total bilirubin 


Rho 


-.077 


-.065 




Sig. 


p = .452 


p = .528 



a Spearman's rho, two-tailed significance p < .008. 



Table 4 

Correlations between neonatal risk factors and near-term regional brain white matter microstructure mean diffusivity (MD), corrected for PMA-at-scan, are reported for the genu, 
splenium (Splen), posterior limb of internal capsule (PL1C), anterior limb of internal capsule (ALIC), thalamus (Thai), and globus pallidus, (GloP). Risk factors include gestational age at 
birth (GA), mean and peak serum C-reactive protein (CRP), and mean serum glucose, albumin, and total (T) bilirubin. 







Genu 


Splen 


LAL1C, 


RALIC, 


LPLIC, 


R PL1C, 


LThal, 


RThal, 


L GloP, 








MD 


MD 


MD 


MD 


MD 


MD 


MD 


MD 


MD 


RGloP, MD 


GA (N = 66) 


r 


.116 


-.176 


.160 
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Correlations among neonatal risk factors. Risk factors include gestational age at birth (GA), mean and peak serum C-reactive protein (CRP), and mean serum glucose, albumin, and total 
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MD demonstrated the strongest associations with neonatal risk factors, 
higher thalamus MD was associated with lower GA (left: r= —.322, 
p = .009; right: r= —.381, p= .002), lower mean serum albumin 
(left: r= -.276, p = .029; right: r= -.385, p= .002), and lower 
mean total bilirubin (left: r= -.293, p = .020; right: r= -.337 p= 
.007). 

4. Discussion 

This study examined the relationship between neonatal physiologi- 
cal risk factors and early brain structure that was observed on MRI and 
quantified using DTI. We found that more than one-third of infants had 
brain abnormalities observed on structural MRI. Infants with brain ab- 
normalities on MRI demonstrated a similar degree of prematurity, how- 
ever they had a higher incidence of BPD and sepsis and higher levels of 
serum CRP over the first 2 weeks of life, compared to those with no 
brain abnormalities on MRI. Evidence from DTI indicated that the thala- 
mus WM microstructure was most strongly associated with signal ab- 
normalities on structural MRI and with neonatal risk factors. 

Advances in neonatal medicine have increased survival rates and 
improved outcome among VLBW preterm infants. The increased surviv- 
al of VLBW preterm infants is associated with increased rate of 
neurodevelopmental impairment, which is currently 3-4 times that of 
the general population (Williams et al., 2010; Spittle et al., 2011). The 
incidence of CP in VLBW preterm infants is estimated to be 10-15% 
(Litt et al., 2005; Williams et al., 2010; Rose et al., 2009; Woodward 
et al., 2012), and is defined as "a group of disorders affecting the devel- 
opment of movement and posture, attributed to non-progressive dis- 
turbances to the developing fetal or infant brain" (Bax et al., 2005; 
Rosenbaum et al., 2007; CDC, 2013). Developmental coordination disor- 
der may affect another 40% of VLBW preterm children of school age 
(Woodward et al., 2012). Cognitive impairment and/or language delays 
have been reported in 20-40% of school-age children born VLBW 
(Anderson, 2003; Litt et al., 2005; Hack, 2006; Larroque et al., 2008; 
Pritchard et al., 2009; Woodward et al., 2009; Aarnoudse-Moens et al., 
2009; Ritter et al., 2014) and it is estimated that half of these children 
will require special educational assistance, (Wocadlo and Rieger, 
2006; Larroque et al., 2008; Johnson et al., 2009; Delobel-Ayoub et al., 
2009; Iwata et al., 2012). It is not currently possible to predict at near- 
term age which children will experience these physical and cognitive 
impairments, limiting our ability to apply early interventions during 
critical periods of development 

Brain MRI scans are commonly performed at near-term age in very 
preterm infants, and was envisioned as a potential opportunity to pro- 
vide early prognosis in preterm infants. However, prognosis based on 
structural brain MRI has been of limited success in predicting risk for 
neurodevelopmental problems such as motor deficits, cognitive delay, 
and language impairment later in life (Kidokoro et al., 2011; Benini 
et al., 2013). DTI of early WM regions may be a more sensitive biomark- 
er of later neurodevelopment than structural MRI (Arzoumanian et al., 
2003; Rose et al., 2007; Rose et al., 2009; De Brui'ne et al., 2013). Semi- 
automated, atlas-based DTI analysis has potential for clinical implemen- 
tation and therefore was used in the current study. Diffusion was ana- 
lyzed in WM regions which were previously suggested to be reliably 
detected by semi-automated, atlas based analysis of DTI at near-term 
age (Rose et al., 2014) and also mediate motor and cognitive function, 
domains known to be affected by preterm birth. 

4.1. Neonatal risk factors and brain abnormalities observed on structural 
MRI 

Brain abnormalities were identified in 38/102 infants, including 21/ 
60 females (35%) and 17/42 males (40%). No differences in GA or BW 
were found, and while infants with brain abnormalities observed on 
MRI had lower Apgar scores at 5 min after birth, the trend did not 
reach significance. The value of Apgar scores for prediction of outcome 



has been debated, but Lie et al. (2010) found a strong correlation be- 
tween low Apgar scores 5 min after birth and diagnosis of CP before 
age 5 in children with both low and normal BW. The incidence of BPD 
was slightly higher in those with brain abnormalities observed on MRI 
(41%) compared to those without (37%) but did not correlate with 
MRI score or the number of signal abnormalities on MRI. 
Bronchopulmonary dysplasia continues to be one of the main risk fac- 
tors for motor impairments following preterm birth (Schmidt et al., 
2003; Vohr et al., 2005; Karagianni et al., 201 1 ; Schlapbach et al., 201 2). 

Infants with brain abnormalities on MRI had significantly higher 
mean and peak serum CRP levels compared to those with no brain ab- 
normalities on MRI. Furthermore, number of signal abnormalities ob- 
served on MRI correlated to serum CRP levels. CRP is a plasma protein 
involved in inflammatory response via rapid response to cell death or 
injury; simultaneously production of other plasma proteins including 
albumin, decreases, thus, CRP levels have been proposed as a biomarker 
for sepsis and general inflammatory response in neonates (Carlo et al., 
2011; Meem et al., 2011; Fan andYu, 2012). One challenge is determin- 
ing an optimal cutoff value for CRP levels for sensitive and specific iden- 
tification of neonatal sepsis (Abdollahi et al., 2012; Boonkasidecha et al., 
20 1 3 ). Using a combination of biomarkers may hold more diagnostic ac- 
curacy than any one biomarker (Wang et al., 2013). Sepsis is a known 
risk factor for poor neurodevelopmental prognosis and cerebral palsy 
among VLBW preterm infants (Stoll et al., 2002; Schlapbach et al., 
2011; Hentges et al., 2013; Mitha et al., 2013; van Vliet et al., 2013). 
Thus, the observed brain abnormalities may be related to inflammatory 
processes that occurred in the neonatal period. 

Abnormal serum concentration of glucose, hormones, and proteins 
are commonly observed in VLBW preterm infants. Specifically, hypogly- 
cemia was shown to be correlated with increased morbidity and mortal- 
ity (Kao et al., 2006; Hays et al., 2006), as well as poor growth (Ramel 
et al., 2013) in VLBW infants. We did not find significant differences be- 
tween serum levels of glucose, albumin, or bilirubin in those with or 
without brain abnormalities observed on MRI. Bilirubin levels observed 
in our population of neonates were not grossly abnormal, and thus our 
assessment of the role of these risk factors in brain development is 
limited. 

4.2. Neonatal risk factors and WM microstructure on DTI 

42.1. Corpus callosum 

Diffusion imaging scalars were analyzed within association fibers of 
genu and splenium of the ccs, regions that are in preliminary stages of 
myelination at near-term age (Brody et al., 1987; Kinney et al., 1988; 
Provenzale et al., 2012). The posterior aspects (splenium) of the ccs 
were previously found to develop earlier than the anterior (genu) as- 
pects (Vassaretal., 2013; Rose etal., 2014). Previous studies also report- 
ed that WM structure in the near-term CC may have prognostic value 
for future motor function (Counsell et al., 2008; Thompson et al., 
2012). Males were found to have significantly higher genu MD, 
corrected for PMA-at-scan. The limited data on sexual dimorphism of 
neonatal CC development suggests that males may have slightly smaller 
or less developed CC at near-term age compared to females (Rose et al., 
2009; Thompson et al., 2012). Male gender also was found to be associ- 
ated with increased risk of adverse outcomes including CP and cognitive 
deficits following preterm birth (Wood et al., 2005; Vohr et al., 2005; 
Rose et al., 2009; Beaino et al., 2010; Morsing et al., 201 1 ). Further inves- 
tigation is needed to explain the physiological and neurological risk fac- 
tors that increase morbidity among VLBW males. 

4.2.2. Internal capsule 

Diffusion was also analyzed in the projection fibers of the anterior 
and posterior limbs of the internal capsule. The PLIC includes motor fi- 
bers of the cortical spinal tract (CST) and has been found to be relatively 
more developed than the ALIC at near-term age (Vassar et al., 2013; 
Rose et al., 2014). Based on the work of Yakovlev and Lecours (1967) 
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and Brody et al. (1987), regions of the CST are believed to be undergoing 
myelination during the near-term period. The CST descends through the 
posterior third of the PLIC and, consistent with this neuroanatomy, anal- 
yses of PLIC DTI at near-term age suggest that it may hold prognostic 
value for future motor development (Arzoumanian et al., 2003; Rose 
et al., 2007; Rose et al., 2009; De Brui'ne et al., 2013). No further associ- 
ations between neonatal risk factors and PLIC MD were found. This may 
be explained by the developmental trajectory of the PLIC, which is 
thought to be particularly vulnerable at the end of the second trimester, 
when developing vasculature and early myelination processes may in- 
crease vulnerability to WM injury (Hiippi et al., 1998; Johnston and 
Hoon, 2000; Neil et al., 2002; Miller and Ferriero, 2009; Liu et al., 2013). 

4.23. Thalamus and globus pallidus 

We found that thalamus WM microstructure demonstrated the 
strongest associations with neonatal risk factors, including GA, albumin, 
and bilirubin. Among VLBW neonates, GA has been found to be inverse- 
ly correlated with cognitive and motor function later in childhood in 
several large studies (Bhutta et al., 2002; Constantinou et al., 2005; 
Odding et al., 2006; Aarnoudse-Moens et al., 2009; Woodward et al., 
2009; Moore et al., 2012). Vulnerability of thalamus WM structure to 
neonatal risk factors could in part explain the impaired motor and cog- 
nitive function associated with preterm birth. 

Two recent studies by Ball et al. observed decreased thalamus vol- 
ume (Ball et al., 2012b) and impaired thalamocortical connectivity 
(Ball et al., 2012a) in infants born preterm, compared to age-matched 
term infants, suggesting that WM impairment may occur due to pre- 
term birth. The thalamus is undergoing significant development and es- 
tablishing cortical connections during the late second and third 
trimester of gestation (Kostovic and Judas, 2010), and this process 
may be disrupted by preterm birth. The thalamus and basal ganglia 
(consisting of the caudate nucleus, putamen, and globus pallidus) are 
known to mediate motor processes and have been found to be vulnera- 
ble to hypoxic ischemic injury in infants (Cowan et al., 2003; Miller and 
Ferriero, 2009; Ball et al., 2012a; Nossin-Manor et al., 2012). Other in- 
vestigators have also reported reduced volumes of these subcortical re- 
gions in preterm infants, compared to age-matched, term-born controls 
(Boardman et al., 2006; Srinivasan et al., 2007). 

We found that higher thalamus MD was associated with lower mean 
albumin and with lower mean total bilirubin. Peak serum bilirubin dur- 
ing the first 2 weeks of life is known to be associated with increased risk 
of neurodevelopmental impairments (Oh et al., 2003; Oh et al., 2010). 
However, aggressive treatment with phototherapy may cause oxidative 
injury to cell membranes and increased mortality (Morris et al., 2008). 
Thus, the effect of high levels of bilirubin and aggressive bilirubin pho- 
totherapy on overall outcome among VLBW infants needs further inves- 
tigation (Tyson et al., 2012), and is difficult to study independently 
(Hintz et al., 2011). 

Recent studies suggest that unbound bilirubin is more informative 
than total serum bilirubin in prediction of neurodevelopmental trajec- 
tories in VLBW infants (Calligaris et al., 2007; Oh et al., 2010). The 
ratio of bound/unbound bilirubin is moderated by the availability of al- 
bumin, as bilirubin binds to albumin, preventing it from crossing the 
blood-brain barrier (Watchko, 2006; Bender et al., 2007). Thus, low al- 
bumin levels in the preterm infant may also be involved in the cytotox- 
icity of bilirubin and may hold prognostic value. Because unbound 
bilirubin is difficult to measure compared to total bilirubin, the ratio of 
total bilirubin/albumin may allow for estimation of exposure of the neo- 
natal brain to unbound bilirubin (Sato et al., 2012). 

The strong association found between albumin and total bilirubin 
levels in the current study suggests that low albumin levels may be pri- 
marily responsible for both increased thalamus MD and lower bilirubin 
levels. Phototherapy to reduce bilirubin levels is administered based on 
BW, and thus infants with lower BW will have lower bilirubin levels 
(Table 5). This may also contribute to the relationship between low bil- 
irubin and high thalamus MD (Table 4). However, further study is 



required to examine the interactions between bilirubin, albumin, and 
sepsis and how they individually contribute to increased risk for 
neurodevelopment impairments. Furthermore, albumin levels are de- 
creased as a result of inflammation, and ischemia-modified albumin 
levels may be used as prognostic tools for sepsis in neonates 
(Yerlikaya et al., 2014; Yakut et al., 2014). 

5. Conclusion 

This study examined the relationship between neonatal physiologi- 
cal risk factors and early brain structure in VLBW preterm infants. DTI 
was assessed to examine evidence of selective vulnerability of brain re- 
gions, which may ultimately serve to predict outcome and guide neuro- 
protective treatment for preterm infants. We used semi-automated 
atlas-based analysis of DTI because of its potential for clinical neuroim- 
aging for neonates. We observed that the thalamus demonstrated the 
most associations with neonatal risk factors, and although these were 
low to moderate correlations, the findings may reflect selective vulner- 
ability of the thalamus at near-term age. Limitations of this study in- 
clude the cross-sectional design, the difficulty in studying risk factors 
independently, and the low incidence of some risk factors affecting 
brain development (e.g. NEC). Sequential neuroimaging during the 
near-term period would be helpful in revealing trajectories of develop- 
ment and how they relate to neonatal risk factors. Furthermore, correla- 
tion between clinical risk factors and regional brain structure and 
microstructure at near term age does not describe the relationship be- 
tween brain abnormality and neurodevelopmental impairment. Results 
may have been influenced by accuracy of automated ROI generation, 
which is better for larger ROIs. As in any DTI investigation, results may 
also be influenced by differences in SNR between brain regions. Future 
studies using different diffusion-weighted analyses and scan sequences 
should be performed to further investigate WM structure in VLBW neo- 
nates and their relationship with neonatal physiological risk factors. Fu- 
ture studies may determine whether associations between neonatal risk 
factors and regional patterns of microstructural development could pro- 
vide early prognostic biomarkers for future function that may ultimately 
guide early intervention aimed at reducing neurodevelopmental 
disabilities. 
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